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Abstract 

Several new relationships between hypergeometric functions are found by com- 
paring results for Feynman integrals calculated using different methods. A new 
expression for the one-loop propagator-type integral with arbitrary masses and ar- 
bitrary powers of propagators is derived in terms of only one Appell hypergeometric 
function F\ . From the comparison of this expression with a previously known one, 
a new relation between the Appell functions F\ and F4 is found. By comparing 
this new expression for the case of equal masses with another known result, a new 
formula for reducing the F\ function with particular arguments to the hypergeomet- 
ric function 3F2 is derived. By comparing results for a particular one-loop vertex 
integral obtained using different methods, a new relationship between F\ functions 
corresponding to a quadratic transformation of the arguments is established. An- 
other reduction formula for the F\ function is found by analysing the imaginary 
part of the two-loop self-energy integral on the cut. An explicit formula relating the 
F\ function and the Gaussian hypergeometric function 2-F1 whose argument is the 
ratio of polynomials of degree six is presented. 
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1 Introduction 



Radiative corrections to different physical quantities needed for the comparison of theo- 
retical predictions with experimental data to be collected with the CERN Large Hadron 
Collider (LHC) and, in future, with an International Linear Collider (ILC) and other col- 
liders are expressed in terms of complicated Feynman integrals. In many cases, radiative 
corrections must be evaluated analytically to achieve reliable accuracies in the calcula- 
tions. The difficulties in calculating Feynman integrals are usually related to the fact that 
they depend on several kinematical scales, i.e. they are functions of several variables. 

Nowadays, one of the most frequently used methods for calculating Feyman integrals 
is based on the Mellin-Barnes integral representation [L2.3J. In many cases, however, this 
method leads to complicated expressions in terms of hypergeometric functions with many 
variables. In order to calculate analytically integrals with several kinematical variables 
and masses, new effective methods are to be developed. Rather promising methods for 
analytic calculations of Feynman integrals are based on recurrence relations. These can 
be recurrence relations with respect to the exponent of a propagator of the integral [I] or 
the parameter of the space-time dimension [5,6, 7J. As was already observed in the one- 
loop case, the solutions of dimensional recurrences are combinations of hypergeometric 
functions [5TIB"f5] . This is also true at the two- loop level [7J. 

As was realized many years ago in Ref. [9], Feynman integrals are generalized hyperge- 
ometric functions. This conjecture was confirmed through the evaluations of specific Feyn- 
man integrals. Some results for Feynman integrals expressed in terms of hypergeometric 
functions may be found in Refs. [Tll2T3lfT0lfTM2l^^ These 
results were obtained using rather different methods, e.g. by directly evaluating the inte- 
grals from their Feynman parameter representations, by applying Mellin-Barnes integral 
representations, by solving recurrence relations, by making use of the negative-dimension 
approach [25], or by using spectral representations. 

As a method for finding relations between hypergeometric functions, the authors of 
Ref. [26] advocated the evaluation of integrals reducible to hypergeometric functions by 
several different methods and the comparison of the results thus obtained. In this respect, 
the evaluation of Feynman integrals may be considered as a rich source for finding rela- 
tions between hypergeometric functions. New transformation and reduction formulae for 
hypergeometric functions were derived by calculating Feynman integrals already a long 
time ago [12]. Several new reduction relations for the Appell hypergeometric functions 
Fi and F 4 obtained by comparing different results for the same Feynman integral were 
presented in Ref. [27] . 

The analytic evaluation of Feynman integrals offers us a unique possibility to find re- 
lations between hypergeometric functions which can be useful in many other applications, 
far away from high-energy physics. On the other hand, the problems emerging when 
evaluating Feynman integrals may become interesting for mathematicians, and their par- 
ticipation in the solution of these problems may lead to essential progress in the evaluation 
of Feynman integrals. 

Our paper organized as follows. In section |2j we present a new result for a one-loop 
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propagator-type integral with arbitrary exponents of propagators and arbitrary masses. 
In section [3j a new formula for the reduction of the Appell function F 4 to the function Fi 
is derived. Setting the masses in the result derived in section 2 to be equal and comparing 
the outcome with a known result, a new formula for the reduction of the Appell function 
F\ to the hypergeometric function 3F2 is obtained. In section HI from the results for the 
one-loop vertex-type integral, a quadratic transformation formula for the Appell function 
Fi is derived. In section |5l from the comparison of results for the imaginary part of a 
two-loop self-energy integral obtained using different methods, a formula for the reduction 
of the Fi function to the Gauss hypergeometric function with a complicated argument is 
obtained. In section El we present a short summary of our results. 



2 New analytic expression for the one-loop propagator- 
type integral 




Figure 1: Feynman diag ram corresponding to the integral Iu^ U2 - 

In this section, we consider the evaluation of the one-loop propagator type integral 
with arbitrary masses and arbitrary powers of propagators, 

W m l> TO 2! s i2) = J md/2 ^ _ pi)2 _ _ p2 y _ m 2], 2 • I 1 ) 

Here and below, it is understood that the usual causal prescription of the propagators is 
used, i.e. l/[k 2 — m 2 ] -h- l/[k 2 — m 2 + iO}. The Feynman diagram corresponding to this 
integral is presented in Figure 1. By using the formula 



a ^i^ 2 r(z/i)r(z/ 2 ) i Q [ax + 6(1 - a;)]^+^ 

the product of the propagator factors can be transformed to an integral over Feynman 
parameters. Changing the integration momentum as q = t + p 2 +x(p\ —P2) and applying 
the formula 

J_L \ = ( ir r (^-i) m 

m d l 2 (t 2 — M 2 ) u y ' (M 2 )-lr(z/)' 1 ' 
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we obtain the following representation for the integral of Eq. ([T]): 

fm 2 m 2 - Sl2 ) = ( m^ r ^+^-i) Z 1 ^rx^Hl-xr- 1 

2 ' ; 1 ; r(z/ 1 )r(z/ 2 ) y [si2X 2 + x(m 2_ m 2_ Si2) + m 2 ] , 1+ , 2 -r 

(4) 

Representing the quadratic polynomial in the denominator as 

S12X 2 + x{m\ — m 2 — S12) + m\ = m\{l — Xix){l — x 2 x), (5) 
and then comparing our integral with the integral representation for the Appell function 

„, ,„ \ T(c) f 1 du u a - l (l -uf- a - 1 

the following result follows: 

r(d) (22 (-l)^r(// 1 + Z/ 2 -f) 

W™i>™ 2 ; s 12 ) - + ^^y^-d/2 

r, ( d d \ 

x Ft I z/i, ui + v 2 - -, v x +v 2 - -; v\ + v 2 \ x_, x+ 1 , (7) 

where 

l+x-y± y/A(l,x, y) 
x± = 2 , 

x = — , y=—, (8) 

with 

A(x,y,2) = (x - y - zf - Ayz. (9) 

For the particular case v\ — v 2 = 1, an expression for the propagator- type integral in 
terms of the Appell function Fx was given in Refs. [TT1I27II29] . 

The integral /i^ 2 (m 2 ,m 2 ; s\ 2 ) is symmetric with respect to the change ^i,m 2 -H- 
u 2 ,m 2 . To understand this symmetry, we first observe that, under the change m 2 -H- m 2 , 
the arguments of the Appell function F± in Eq. (J7J) transform as 

x_ -> X ~ , x+ ->■ X+ , (10) 
x_ — 1 x + — 1 

and then, applying the formula (see, for example, Refs. |28f30f31j ) 

F x (a, fi, /?', 7, w, z) = (l- w)-P(l - z)^F x (7 - a, /3, /?', 7, -^r) , (H) 

w — \ z — l 

we return to the initial expression on the right-hand side of Eq. (J7J) . 
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3 Relations between the F\ function and other hy- 
pergeometric functions 

In Ref. PQ, by exploiting the Mellin-Barnes integral representation, the following analytic 
expression for the considered integral was derived: 



( _1 

d d d S\2 mf 

2^' 2 
2 m 2 



+ '4J -i^r^^?2-2-^ +1 '4 i 4 l> - (12) 

A hypergeometric representation in terms of Lauricella functions for the one-loop integrals 
corresponding to diagrams with an arbitrary number of external legs was presented in 
Refs. [2113] . 

Comparing Eqs. (J7J) and ({12]) . we arrive at the following relation: 

n f d d 

Fi I z/i, Ui + u 2 - -, vx + v 2 - -; V\ + V 2 \ x H 

r (I - i/Q rfo + v 2 ) ( dd d 

«, YU-f) I>! + 1/ 2 ) ^ / d d d \ 

r(i/i)r (i/i + 1/ 2 - f ) v 2 2 2 / 

Here, a;± are given by Eq. (JSJ). With the help of the relation given in Ref. [32J (see p. 102), 
Fa ( a,P;7,P\ r|- r, r|- ) = [(l-x)(l-y)] oc F 1 (a,'y-P, l+a-7, 7; x, xy), 



(l-x)(l-yY (l-x)(l-y) 

(14) 

the second Appell functions F4 on the right-hand side of Eq. ( 1T31) may be expressed in 
terms of the Appell function Fx. Therefore, the following relation holds: 

F 4 (a, (3, (3', a — (3' + 1; x, y) 

r(^)r(/3-a + /30r(a-^) . a _ ^ 
r(«)r(/3)r(/3'-«) y K+ ' 

xF 1 (p-a + P',l-a,P,P'-a + l; ffiZifh £Zfh ) _ ( 15 ) 
V x x — x. J 
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In Ref. [3], an expression for the integral Iv\v2 (mf, ml; s 12 ) in terms of the Kampe de 
Feriet function was derived: 



j(d) 



^' U2 {ml,ml; s 12 



A r l;0;0 



) = (-ir +1 *(m|)*- 



r(z/i + 1/ 2 ) 

(^ + ^-d/2:l,l),(^ : 1,1) : (y 2 : 1) 
(i/i + v 2 ■ 2, 1) 



^12 _ ■ j_ 

2 ' 2 
m 2 m 2 J 



(16) 



Comparing this relation with Eq. (J7|), we obtain the following reduction formula: 



2:1:0 



(a:l,l),(^:l,l):(^:l) 
(i/i + i/ 2 : 2, 1) 



Fi(vi_,a,a,vi + v 2 \ Z-,z+) 



(17) 



where 



z± 



x + y ± i/(x + y) 2 — Ax 



For the case of equal masses ml 
Ref.|I]: 

I^Jm 2 ,m 2 ; s 12 



m\ 



2 v 7 

m 2 , the following expression was derived in 



^1 1^2 s 



x — rv — ; — ^ — 3 ^ 2 
r(^i + z/ 2 ) 



I/l, I/ 2 , ^1 + ^2 



., "Zj L i z 2 , 
2 ' 2 ' 



Am 2 



Comparing this formula with Eq. ([7j) taken at m\ = ml = m 2 , we obtain: 

= 3-^2 



F\ ( a, /3, /3; 7; x — V x 2 — 2x, x + y 



which may be rewritten as: 

F 1 [ a,/3,/3;7; x, 



x 



X 



3-^2 



2x 



a, 7 — a, j3 ; 

7 7+1 



a, 7 — a, /3 ; x 
2 



, / H 1 

7 7+1 . 
2' 2 ' 



X 



2> 2 ' ^(x 



(19) 



(20) 



(21) 



We verified numerically the correctness of this relation setting a = 1/2, /3 = 2, 7 = 3/2, 
and a; = —1/2, and keeping 600 valid digits in the calculations performed using the 
computer algebra system Maple. 

To the best of our knowledge, there is no such a relation in the mathematical literature, 
i.e. Eq. ( l2~Tj) extends the number of known reduction formulas for the Appell function 
In Ref. [33] , a relation between Fi with the same arguments and the Gauss hypergeometric 
function 2 Fi is given. That relation corresponds to a particular case of our Eq. (|2~Tj) . taken 
at 7 = 2a. We would like to recall that the only known transformation of the Appell 
function F\ to the hypergeometric function %F 2 was known for the case when x = —y. 

For the Appell hypergeometric function F\, the following relation holds: 



Fi(a,/3,/3;7;x,-x) 



3-^2 



a q+1 a . 
2 ' 2 ' P ' ™2 

7 7+I . •* J 

2' "2" ' 



:22i 
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In the case when su 
rid) 



m,, the relation holds. Therefore, the integral 



(mf, m 2 ,; S12) may be expressed in terms of the hypergeometric function 3 F 2 as 



si2=mf — rn\ 



_ (-1)^1+^(1/! + l/ 2 -f) 



3-^2 



T' _ 2~ ' 1 2 ~~ 2 ' 

1^1+^2 ^1+^2 + 1 . 
2 ' 2 ' 



1 _ !^ 

2 

m 2. 



(23) 



This formula demonstrates that simplifications of Feynman integrals may also take place 
for specific values of masses or momenta that are more general than just zero or on-shell. 



4 Quadratic transformation for the Appell function 




q-pi 



q-P2 



q-P3 



<d) 



Figure 2: Feynman diagram corresponding to the integral I 

In this section, we find relations for the Appell function Fi by comparing the results 
of different calculations of the one-loop vertex-type integral 



r(d) / 2 2 2 \ 

P 3 '{m^m^m^p^pjupjk) 



1 



i^d/2 — pj-) 2 — m^] [(g - Pk) 2 -rn 2 k ][(q -pi) 2 -m 2 y 



(24) 

corresponding to the Feynman diagram shown in Figure 2. We consider this integral with 
a particular set of arguments, namely 1^(0, m 2 , 0; 0, S13, S12). In the case when s 12 < m 2 
and S13 < m 2 , from its representation as an integral over Feynman parameters, we derive 
the following analytic expression: 



4 d) (0,m 2 ,0; 0,s 13 ,s 12 ) 



1 r i 



—dx\d%2 r (3 — |) X\ 



1 

m 



2 /Jf(0,m 2 ; 0) F l ( 1,1,2 



Jo [s 13 Xi - S13 + X 2 {m 2 + S13 - S12 - X1S13 + X1S12)] 3 2 

d cf S12 — S13 S12 
2' 2 : 



rrr 



A ( f(0,0; s 13 ) 



^1 



1, -3— , -^12 ~ g 13 



tZ-2; 



(25) 
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where F\ is the Appell hypergeometric function [28] defined by the series 

P/ x {a) k+ i{b) k (b')iw k z l 

Fx (a, b,b;c;w,z) = > — — , 26 

k j^ (c)k+l kill 

and (a)k = T(a + k)/T(a) is the so-called Pochhammer symbol. In our case, the Appell 
function has a rather simple integral representation, viz. 

Therefore, by using Eq. (I25I) . one may obtain a result for the integral 1% in terms of the 
Appell function F\. The result in terms of function F\, was previously obtained in Ref. [8] 
and later on in Ref. (21]. In d = 4 space-time dimensions, the result for the integral 1^ 



in terms of the function F 3 was given in Ref. [20J. 

The integral 1^(0, m 2 , 0; 0, s i3 , Si 2 ) may also be evaluated by another method, based 
on difference equations with respect to the space-time dimension d. The method of de- 
riving dimensional recurrences is described in detail in Refs. 06]. In the case under 
consideration here, we have 

T (d+2) /n 2 n n \ 2m2s 13(si2 ~ m 2 - S 13 ) ( d ) . 2 

'{0,771 ,0; 0,si 3 ,si 2 ) = — -, r^-. — -t — I3 '{0,m ,0; 0, si 3 , s 12 ) 

{S12 — S13) [a — i) 

m 2 {d) 2 _ m , s 13 {s 12 - gig ~ 2m 2 ) (d) 

+ (m2gl2 / m2gl \t;; 2gl 9 3 "^ Iff (0,m 2 ;, 12 ). (28) 
(S12 - si 3 y{d - 2) 

Denoting its non-homogeneous part as R^ d \ we rewrite Eq. (I28I) as: 

/r 2) (0,m 2 ,0; 0, s 13 , s 12 ) = 2r f Sl3(g12 T^ 3 ~ f } 4 rf) (0, m 2 , 0; 0, s 13 , s 12 ) + flW. (29) 

(sia - si 3 ) 2 (a- 2) 

The solution of this equation reads 



/f(0,m 2 ,0; 0, Sl3 , Sl2 ) = -42YC £ ( Sl2 , Sl3 )-^-^ f;%i2(^») 



r(^) £K iZ ' ^ 2a 

where 



fc=0 



(30) 



m 2 si 3 (si 2 - si 3 - m 2 ) 
(T= 7 • (31) 
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An arbitrary periodic function C e (si 2 , $13) emerging in the solution may be found from 
the following differential equation with respect to the variable si 2 : 







4 d) (0,m\0;0,s 13j s 



ds 



12 



12 



{d - 2)(sia - S12 + 2m 2 ) - 2m 2 
2(m 2 + si 3 - si 2 )(si 3 - S12) 



4 d) (0,m 2 ,0;0,s 13 ,s 12 ) 



+ 



+ 



(rf-3)(m 2 + s 13 -2s 12 ) 
(m 2 - si 2 )(m 2 + si 3 - s 12 )(s 12 - s 13 ) 



lif(0,m 2 ;s 12 ) 



(d-2) 



/ff(0,m a ;0) 



2(m 2 - si 2 )(m 2 + S13 - s i2 

(d-3) ,,,/, 



J n (0, 0; s 13 ). 

(m 2 + S13 - S 12 ){s 12 - S13J 

Substituting Eq. f l3"Uj) into Eq. f[3"2j) . we arrive at the following equation: 







■C e (s 12 , si 3 ) + 



S12 - S13 - 3m 2 



-C £ (s 12 ,s 13 ) 



0. 



(32) 



(33) 



aSi 2 ^ tV ""'" ±0 ' ' (s 12 - S13 - m 2 )(si 3 - S12) 

Taking into account the boundary condition of the integral at s\ 2 = 0, the solution of this 
equation is 

C £ (SX2, S13) = 0. 

Substituting explicit expressions for the propagator integrals I 2 into Eq. (130]) leads to 
the following expression: 

lf(0,m 2 ,0; 0,s 13 ,s 12 ) = 



(m 2 - si 2 )si2 + (^ 2 + s 12 )s 



2si 3 (s 13 - S12 + m 2 )(s 12 - m 2 ) 

S13 - S12 r (d) /n _ 2 



^/;f(0,m 2 ; 0)F 1 [ - 1 ^,-,1,~, 



d-2 1 
2 ' 2' 



Am S\ 2 m 
2 ' (S12 — m 2 ) 2 ' 



+ 



+ 



2(si3 - s 12 + m 2 )s i3 

(gig ~ -Si3 ~ 2m 2 ) 
2m 2 (si 3 - s 12 + m 2 ) 



I7(0,m 2 ; 0) 2 F X 
4? (0,0; s 13 ) 2 Fi 



l,^;-m 2 
<r. 



2 ■ 



1 rf - 2 ■ c 
1, — , s 13 

d-1 . ~ A 

— ; 4a 



Comparison of the obtained result with Eq. (I25I) leads to the relation 



d d 

Ft ( 1,1,2--, ~;u,z 



00 



+ 



(u + zu — 2z) 



2(u-z)(l-u) 
(d-2 1 d 

Fx (— — ,1 - -; 



2^1 



1 ^2 . 

2 , 



§; (w-«)(w-l)_ 



-4z 



a 



(34) 



2(w-z)(l-w)(l-^) * A V 2 '~'2'2' (w-2)(w-l)' (*-l)V ' (35) 

The arguments of the function F\ on the right-hand side of Eq. (I35I) are connected with 
the arguments of the F\ function on the left-hand side by a quadratic transformation. 
Therefore, Eq. (135|) is the analogue of the quadratic relation for the Gauss hypergeometric 
function 2 F\. To the best of our knowledge, there is no such a relation in mathematical 
literature. 
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5 New relation between Appel function F\ and hy- 
pergeometric function 2^1 




Figure 3: Feynman diagram corresponding to the integral Jg . 

In this section, we find a relation between the F\ and 2 -Pi functions by comparing the 
results evaluated by two different methods for the imaginary part of the integral 

{d) = r r d%d% 1 

~JJ {m d ' 2 ) 2 {k 2 - m 2 )^ - k 2 ) 2 - m 2 ){{k 2 - q) 2 - m 2 )' 1 ' 

The Feynman diagram corresponding to this integral is presented in Figure 3. In Ref. [TJ, 



the difference equation with respect to d for the master integral was derived. From 
this equation, the difference equation for the imaginary part of may be obtained. It 
reads: 

12x 2 {d + l)(d - l)(3d + 4) (3d + 2) Imjf +4) 



-4m 4 (i 2 - 3)(x 4 - 42x 2 + 9)(d - l)d ImJ. 



-4m 8 (x 1 



D 2 (x 2 



(d+2) 
3 



9) 2 ImJ^ = 



(37) 



where a? = q/m. The solution of this equation for the imaginary part of Jg was presented 
in Ref. [7J and reads: 



ImJ. 



(d) 



-4 vr 2 V3 m M - 6 
T (d - 1) [x 2 + 3) 



(x 2 - 9) 2 
27 



2^1 



l,l; xV-9) 2 
6 d. 

2 ' 



(x 2 + 3f 



(3f 



An analytic expression for ImJg may also be obtained by using another method. In 
Ref. [34J, a one- fold integral representation for the imaginary part of the two- loop sun- 
rise integral with arbitrary masses was derived. For our case, where the masses of all 
propagators are the same, the imaginary part on the cut is given in Ref. [31] and reads: 



ImJ. 



(d) 



-7T 



r 2 (¥) 



(g2)|-ir 2 (d-2) J 4m 2 



(q—m) 2 



d9 



-(A(6,m 2 ,m 2 )A(6,q 2 ,m 2 )) 



d-3 

2 



(39) 
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where A(x, y,z) is defined in Eq. (jUJ). Changing the integration variable in Eq. ([59]) as 
9 = 4m 2 + (q — 3m) (g + m)/3 leads to the following expresssion: 



ImJ. 



('J) 



-* T 2 (¥) (ff-3m)(g + 



m 



( g 2)|-ir 2 (d-2) 2m 



[(g — m)(g + 3m)(g + m) 2 (g — 3m) 2 ] 
1 



d-3 
2 



d-3 
2 



[1 



(q+m)(g— 3m) 
4m 2 



(40) 



As follows from Eq. ([6]), the integral on the right-hand side of Eq. (140]) is proportional to 
the Appell function iq. Using Eq. (jfJJ) and comparing Eqs. ( l4"0j) and fl38|) . we arrive at the 
following relation: 



tZ-13-dl (x + l)(x-3) (x + l)(x-3) 

"l - '^ - '2' ; (z-l)(x + 3)' 4 



2^ 

(x 2 + 3) 



16x 2 (x + 3) 
27 (x + l) 2 



21 



[(x + 3)(x-l)] : 



2^1 



2 . ™2 / ,y»2 
3> 3 > x l x 
d . 
2 ' 



9) 



21 



x 2 + 3) J 



(41) 



It is interesting to note that, at d = 2, the Appell function iq in Eq. fl4"T]) may be 
expressed in terms of the 2 F X function with the help of the equation (see, for example, 
Ref. EH) 



Fx{a,b,b',b + H,w,z) = (1 -z)- a 2 F 1 
This leads to the relation: 

^-3)(x + l) 3 ] _ v/3(x + 3)(x - l) 3 



a, b ; w — z 



b + b'-i 



iFx 



i i . 

2' 2 ' 

1; (:r + 3)(:r-l) 3 



{x 2 + 3) 



2^1 



i,f; xV-9) 2 
1 ; (x 2 + 3) 3 



(42) 



(43) 



The hypergeometric function 2 Fi on the left-hand side of this equation is proportional 
to the complete elliptic integral of the first kind. Relations between hypergeometric 
functions with parameters 1/2,1/2,1 and 1/3,2/3,1 but with arguments different from 
that in Eq. (|43p were first derived by Ramanujan in Ref. [35] . 



6 Conclusions 

In this section, we briefly summarize the most important results obtained in this paper 
and point out some topics which may be of interest for future investigations. Specifically, 

• a new analytic expression for the one-loop propagator-type Feynnman integral was 
derived; 

• a new formula transforming the Appell function F\ to the hypergeometric function 
3 F 2 was presented; 
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• a new formula transforming the Appell function F4 to a combination of two F\ 
functions was found; 

• a new formula connecting the Kampe de Feriet function and the Appell function Fi 
was obtained; 

• a formula for the quadratic transformation of the Appell function F\ was given; 

• a new relation between the Appell function F\ and the elliptic-type Gaussian hy- 
pergeometric function 2 F\ was established. 

We expect that the formulae found for the abovementioned hypergeometric functions 
will be useful, in particular, for finding relations between transcendental numbers. Such 
relations may be found, for example, via the expansions in e — (4 — d)/2 of the presented 
results. These relations may also be useful for simplifying the e expansions of hypergeo- 
metric functions involved because some of them may have simpler integral representations 
than the others. 
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